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Abstract 

The exponential cubic B-spline functions together with Crank Nicolson are 
used to solve numerically the nonlinear coupled Burgers’ equation using collocation 
method. This method has been tested by three different problems. The proposed 
scheme is compared with some existing methods. We have noticed that proposed 
scheme produced a highly accurate results. 


1 Introduction 


The purpose of this paper is to apply the exponential B-spline collocation method 
to the coupled Burgers equation system. The Coupled Burger equation system in the 
following form 


dU d 2 U . dU . /rrTM 

——- + k\U ——b k-2(UV) x — 0 
ox z ox 

d T V dV 

7TT + k ' V lT + H UV ) X = 0 
Ox z ox 


dt 

dV 

~dt 


( 1 ) 


where k \, k 2 and k% are reel constants and subscripts x and t denote differentiation, x 
distance and t time, is considered. Boundary conditions 


U(a,t) = fi(a,t), U(b,t ) = f 2 {b,t) 

V (a, t) = gi(a,t), V(b,t ) = g 2 (b,t), t > 0 
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and initial conditions 


(3) 


U(x, 0) = f(x) 

V(x, 0) = g(x), a < x < b 

will be decided in the later sections according to test problem. 

Various methods are used solve the nonlinear coupled Burgers’ equation numerically; 
which is suggested by [1] firstly. Fourth order accurate compact ADI scheme [2], A 
chebyshev spectral collocation method [3], A meshfree technique [3], the Fourier pseu- 
dospectral method [5], the generalized two-dimensional differential transform method [6], 
cubic B-spline collocation method [7], generalized differential quadrature method [SJ, a 
robust technique for solving optimal control of coupled Burgers’ equations [D], a differ¬ 
ential quadrature method [TO] . Galerkin quadratic B-spline finite element method [TT] . 
a fully implicit hnite-difference method [T2] , a composite numerical scheme based on fi¬ 
nite difference [13], an implicit logarithmic finite difference method [TT] . modified cubic 
B-spline collocation method [T5] was applied to obtain numerical solution of nonlinear 
Coupled Burgers system. There are not many articles about exponential cubic B-spline 
method for the solving nonlinear differential equation system. 

The exponential splines and exponential B-splines are defined as a generalization of the 
well-known splines and B-splines by McCartin |16l fT71 [IS ]. He has also showed a reliable 
algorithm by using the exponential spline functions to solve the hyperbolic conservation 
laws McCartin[19j. The exponential B-splines include a free parameter which cause to 
have different bell like piece wise polynomial. The best free parameter is determined 
for the exponential B-spline functions for solving the differential equations. The use of 
the exponential B-splines is not as common as the well known B-splines. There are a 
few studies existing to use exponential B-splines for build up numerical methods. The 
singularly perturbed boundary value problem has been solved based on the collocation 
methods with the exponential B-splincs [2D, 2T[[22J. Very recently, Exponential B-spline 
collocation method is applied to compute numerical solution of the convection diffusion 
equation [23]. 

The paper is organized as follows. In Section 2, some details about exponential cubic 
B-spline collocation method are provided. In Section 3, the initial states are documented. 
In section 4 , numerical results for three different problems and some related figures are 
given in order to show the efficiency as well as the accuracy of the proposed method. 
Finally, conclusions are followed in Section 5. 

2 Exponential Cubic B-spline Collocation Method 

Let 7r be partition of the problem domain [a, b] defined at the knots 

7T : a = x 0 < Xi < ... < xn = b 

with mesh spacing h = (b — a)/N. The exponential B-splincs, Bi(x), with knots at the 
points of tt can be defined as 
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Bi(x) = < 


where 


(xi -2 - x) - - (sinh(p ( Xi -2 ~ x))) ] [Xj_ 2 , 24 -iJ 

ai + h (xi -x)+c 1 exp (p (xi - x)) + d 1 exp (-p (x» - x)) [xj_i, x*], 
oi + h(x- + ci exp {p (x - x^) + d 1 exp {-p (x - x t )) [x h x i+1 ], 

(x - x i+2 ) -(sinli (p (x - x i+2 ))) 

p 


b 2 

0 


[•£*+1) *^*+2] j 

otherwise. 


(4) 


a\ = 


Cl = 


di = 


phc , p 

? 61 = 2 


phc 


c(c - 1) + g 2 
(phc — s)(l — c) 


b 2 = 


p 


2 (phc 


1 f exp(—ph)(l — c) + s(exp(—ph) — 1) 


4 V (phc — s)(l — c) 

1 h exp (ph)(c — 1) + s(exp(ph) — 1) 


(phc — s)(l — c) 

and s = sinli (ph), c = cosh (ph), p is a free parameter. When p — 1, graph of the 
exponential cubic B-splincs over the interval [0.1] is depicted in Fig. 1. 


Fig.l: Exponential cubic B-splines for p = 1 over the interval [0.1] 

An additional knots outside the problem domain, positioned at x_i < x$ and Xn < 
x tv_)_i are necessary to define all exponential splines. So that {£>_i(x), B 0 (x), • • • , B N+ i(x)} 
forms a basis for the functions defined over the interval. Each sub interval [xj,Xj + 1 ] is 
covered by four consecutive exponential B-splines. The exponential B-splines and its first 
and second derivatives vanish outside its support interval [xi^ 2 ,x i+2 \. Each basis function 
Bi(x) is twice continuously differentiable. The values of Bi(x), B { (x) and B i (x) at the 
knots Xi can be computed from Eq.Q are shown Table 1. 
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Table 1: Values of Bi{x) and its principle two 
derivatives at the knot points 


X 

%i —2 %i— 1 %i +1 %i +2 

Bi 

b: 

Y s-ph “ s-ph, ~ 

2 {phc — s) 2(phc — s) 

0 P(l-c) Q p(c- 1) Q 

2 (phc — s ) 2 (phc — s) 

2 2 2 
ps ps ps ^ 

2 (phc — s) phc — s 2 (phc — s) 


An approximate solution Un(x, t) and VN(x,t) to the analytical solution U(x,t ) and 
V (. x, t) can be assumed of the forms 


N+l N +1 

U N (x,t ) = Y SiBi(x), V N (x,t) = Y ^i B %{x) ( 5 ) 

i=—1 i=—1 

where <5* are time dependent parameters to be determined from the collocation method. 
The first and second derivatives also can be defined by 


N+l N+l 

U' N {x,t) = Y 5 i B i(x), V' N {x,t) = Y <t>tBi{x) ( 6 ) 

i =—1 i =—1 

N+l N+l 

u'+x,t) = y^ b '^ x )^ v nM-Y^ b ^ 

i =—1 i =—1 

Using the Eq. (jSJ), (jBJ) and Table 1, we see that the nodal values U t , Vi, their Erst 
derivatives U[, V( and second derivatives U", V" at the knots are given in terms of pa¬ 
rameters by the following relations 


Ui 

U[ 

U? 

V- - 

VI- 

v!' 


tt! \ s—ph s—ph 

U(X,A = 2 (phc-V- 1 +k+ 2(phc - s) S<+u 
TTU Pil-C) X , P(C- 1 ) x 

U {x, ' t] - 2^—V'- 1 + 2 (phc-8+' 

,w (Xi , t) = 7Tr l±^ Si _^+^ Sl+ 


2 (phc—s) phc — s 2 (phc- 

YM = £++*»+ P(c_1) ' 


2 (phc — s ) 


2 (phc — s ) 


J i +1 


T r„( + 1 P 2S / P 2S . . P 2S 

v (■ Xi,t ) = - - - 7 - + 


2 (phc — s ) 


phc — s 1 2 (phc 




(7) 


The Crank-Nicolson scheme is used to discretize time variables of the unknown U 
and V in the Coupled Burger equation system which is given (jT]) , we obtain the time 
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( 8 ) 


discretized form of the equation as 


U n+1 - U n U2+ 1 + u x 


At , 2 

V n+r_ yn yn+1 + y 


+ fcl W’ +l + ( UU : 


+ A); 


(f/^r 1 + (UV)l 


+ , Avv x r +l2 +(W x n r + h Av)A l2 + ( UV)\ 


= 0 


= 0 


At 2 2 ' 2 

where f/ n+1 = U(x,tn + At) and V n+l = V(x,tn + At). The nonlinear term (UU x ) n+1 , 
(V\ 4) n+1 and (UV)™ +1 in Eq. (jHJ) is linearized by using the following form 


(UU x ) n+1 = U n+1 U% + U n U2 +1 ~ U n U% 

\w x ) n+1 = V n+1 V™ + V n V™ +1 - V n V£ 

(UV)A 1 = ( U x v ) n+1 + ( uv x ) n+1 

= UA l V n + U™V n+1 - U£V n + U n+1 V£ + u n v™ +1 - U n V? 


(9) 


Substitution the approximate solution (J7J) into (JHJ) and evaluating the resulting equa¬ 
tions at the knots yields the system of the fully-discretized equations 


^ml5^\ + ^ m 20m^ + ^m3^m +i + ^m40m +i + ^m5^"+l+ Z/ m6 ( / ) m"+ L l — u m7^m-l^ u rnS^m~^ u rn9^m+l 

( 10 ) 

and 


m+1 


rn+1 


m+1 


cn+1 


m+1 


+ + +^ml30m + ^ml4<5m+l + * / ml50m+l — v m70m-l + z/ mS ( Kn Jrl ' m90m+l 

(U) 

where 

z/ m i = ^ — + hK 2 + A; 2 A 2 ^ ay + (Ay Ad + k 2 L\) ( /3 1 — y, 

^m 2 = {k\K 2 ) ay + (k 2 Ki) A 

Vms — d- k\K 2 + k 2 L 2 S j a 2 — y 2 

^m4 = (AyAT 2 ) a 2 

n m5 = + Ay A" 2 + k 2 L^j ay - (Ay Ad + A; 2 Ai) - y x 

^m6 = (hK 2 ) ay - (k 2 K{) A 


Vm7 = At a i + y x 
Z/m8 = Xt a 2 7 2 
,y m9 = + 7l 

^mio — (k 3 L 2 ) ay + (k 3 L 3 ) / 3 1 

^mii = (At d" k\L 2 + k 3 R 2 ) ay + (AyZa + k 3 R 3 ) / 3 1 — y x 
^rn 12 = ( k 3 L 2 ) a 2 

v m i2 — (At d - k\L 2 + k 3 R 2 ) a 2 — y 2 
^m.14 = (k 3 L 2 ) ay — ( k 3 L 3 ) 

Z7ni5 — (At d- k\L 2 + k 3 K 2 ) ay — (AnAi + k 3 K\) /3 X — y x 
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Ai — (XlS^-l + «2<^m + «3<^m+l -^1 “ a 10m-l + a 24>m + a 30m+l 

k 2 = ^8 n m _ 1 + /3 2 c + ftC + i A 2 = /V™-1 + Pifm + ftC+1 


Oi i — 


s — ph 
2 (phc — s) 


, <22 — 1 , 


o P(!-c) fl P(c-l) 

Pi — o i _ „\ ) P 2 ~ 


2 {phc — s ) 


7i = 


2 

ps 


72 = 


2 (phc — s ) 
2 

p S 


2 {phc — sY ' 1 phc —s' 

The system with (TTOT) and (TTT1) can be converted the following matrices system; 

Ax” +1 = Bx" 

where 


( 12 ) 


A = 


B = 


V ml 

0 


^m3 

^m4 

^m5 

^m6 



u ml2 

^ml3 

^ml4 

Vmlb 



V ml 

^m2 

^m3 

V mA 

^m5 

^m6 

VmlO 

^ m \ 1 

Vml2 

Vml3 


^ml5 



^ml 

v m2 

V m 3 

^m4 ^m5 ^m6 



^mlO 

^mll 

Vml2 

^rn 13 ^ml4 ^ml5 

V m 8 

0 

^m9 

0 




7 0 

V m 8 

0 

Vm9 




^m7 

0 

Vm8 

0 

V m9 

0 


0 

^m7 

0 

^m8 

0 

^m9 




^m7 

0 

^m8 

0 ^m9 0 




0 

^m7 

0 

^m8 0 ^m9 



The system (1T2|) consist of 2 N + 2 linear equation in 2 N + 6 unknown parameters 
x n+1 = (5”!', , <5q +1 , 0q +1 , ..., (j>nX i) )• To obtain a unique solution an additional 

four constraints are needed. By imposing the Dirichlet boundary conditions this will lead 
us to the following relations; 


h_i = (/i(a, t) — a 2 5o — a 3 Si)/ai 
0-1 = (gi(a,t) - a 2 0o - «30i)/ai 
8n +1 = (/2(fr, t) — ai^TV-i — a 2 h7v)/<23 
0W+1 = (fi , 2(&, £) — cti 0 JV _ 1 + a 2 (p N )/a 3 
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3 The Initial State 


Initial parameters <5° 1; 0° l5 5°, <$j,..., b° N+1: (j)°N + i can be determined from the initial 
condition and first space derivative of the initial conditions at the boundaries as the 
following: 


s — ph 


U°(a, °) = 1 + + WTZT^Z~ S i 


s - ph 0 


2 (phc — s ) 


2 (phc — s ) 


U°(x m , 0) = + S° m + A. - 


2 (phc 


m +1 ’ 


2 (phc — s 

u«(b, o) = c/» = J~ ph +< + ,.T pft 


m = 1, 2,IV — 1 


(14) 


2 (phc — s ) 


2 (phc 


N+l 


and 


V°(a, 0) = 


ph 


V°(x m ,0) = 
V°(b, 0) = 


2 (phc — s ) 
s — ph 


£, + *8+ s_pft 


2 [phc — s ) 


C-i + 


2 (phc — s) 1 
o , s ~P h I o 


+ 


2 (phc — s) 


-Vrt+lJ 


m = 1, 2,IV — 1 


(15) 


s - ph , n , n s - ph , n 

;0jv-i + 0at + —;T07V+1 


2 (phc — s ) 


2 (phc — s) 


The system (IITD which is constituted for initial conditions consists IV + 1 equations 
and IV + 3 unknown, so we have to eliminate and 8% +1 for solving this system using 
following derivatives conditions 
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o 

-i 



2 (phc — s) 
p(l - c) 




o _ 

N+l ~ 


6 


o 

N-l 


2 (phc — s) 
P( 1 - c) 




and if the equations system is rearranged for the above conditions, then following form is 
obtained 


^ s—ph 

phc—s 

s—ph -i s—ph 

2{phc—s) 2 (phc—s) 


s—ph 
2 (phc—s) 


l 

s—ph 

phc—s 


s—ph 
2 (phc—s) 

l 


<5? 


77' _ TV 

U 0 p(l-c) U 0 




U'n -i 




_s-ph_Tj/ 
p(c— 1) -W 


which can also be solved using a variant of the Thomas algorithm. As the same way, from 
the system ([15]) . and can be eliminated using 




2 (phc — s ) 
p(l - c) 




2 (phc — s) 
p(l - c) 
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conditions and the following three bounded matrix is obtained. 


2 s—ph 

phc—s 

s—ph i s—ph 

2 (phc—s) 2 (phc—s) 


K 


r TV s~ph TV 1 

v 0 p(i-c) 

v; 

s—ph i s—ph 

2 (phc—s) 2 (phc—s) 

s—ph 
phc—s 


< f > %-1 


Vir-i 

TV s~ph TV 

L V N p(c-l) V N J 


4 Numerical Tests 

In this section, we will compare the efficiency and accuracy of suggested method 
problem. The obtained results will compare with [7], |10j . [15] and [25], while p changes. 
The accuracy of the schemes is measured in terms of the following discrete error norm 

Eoo 

L 00 = \U- C/ivL = max \U j - (U N )]\ . 

Problem 1) Consider the Coupled Burgers equation system (|T[) with the following initial 
and boundary conditions 

U(x, 0) = sin(x), V(x, 0) = sin(x) 

and 

U(—7T,t) = U(7T,t ) = V(—7T,t) — V(n,t) — 0 

The exact solution is 

U(x,t ) = V(x,t) = sin(x) 

We compute the numerical solutions using the selected values k\ = —2, k 2 = 1 and 
k ?J — 1 with different values of time step length At.In our first computation, we take 
t = 0.1, At = 0.001 while the number of partition N changes. The corresponding results 
are presented in Table 2 a. In our computation, we compute the maximum absolute 
errors at time level t = 1 for the parameters with different decreasing values of t. The 
corresponding results are reported in Table 2 b. In both computations, the results are 
same for U(x,t ) and V(x,t) because of symmetric initial and boundary conditions. And 
also we correspond the obtained numerical solutions by different settings of parameters, 
specifically for those taken by [15] in Table 2 c for N = 50, At = 0.01 and increasing t. 
And also in Table 2, we present the rate of convergence in space which is clearly of second 
order. 


Table 2 a: 

Lqo Error norms for t = 0.1, A t = 0.001, U(x,t ) = V(x,t) 




Present (p — 1) Present (Various p) 

m, (a = o) 

p5] (Various A) 

N 

= 200 

0.01489 x 10" 5 0.00121 x 10~ 5 

0.74326 x 10~ 5 

0.00079 x 10" 5 



(p = 0.00004330000) 


(A = -1.640 x 10 -4 ) 

N 

= 400 

0.00372 x 10" 5 0.00044 x 10~ 5 

0.18534 x 10" 5 

0.00006 x 10" 5 



(p = 0.00012611302) 


(A = -4.087 x 10" 5 ) 
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Table 2 b: Error norms for t = 1, N = 400, U(x,t ) = V(x,t) 



Present (p = 1) 

Present (Various p ) 

m, (a = o) 


(Various A) 

At = 0.01 

1.8194 x 10“ 5 

0.00247 x 10“ 5 
(p = 0.00010099997) 

1.08691 x 10” 5 

(A 

0.00131 x 10~ 5 
= -5.896 x 10" 5 ) 

At = 0.001 

1.5159 x 10“ 5 

0.00309 x 10“ 5 

(p = 0.00021660000) 

1.10393 x 10” 5 

(A 

0.00036 x 10~ 5 
= -5.992 x 10" 5 ) 


Table 2 c 

: Loo Error norms 

for At = 0.01, N = 

50, U(x,t ) = V(x,t) different t. 




Present (p — 1) 

Present (Various p) 

HD] 

HU 

m 

t = 

0.5 

7.9881 x 10" 4 

3.4770 x 10~ 4 

1.51688 x 10 -4 

2.26627 x 10” 5 

1.103080984 x 10~ 4 

t = 

1.0 

9.6837 x 10" 4 

4.2166 x 10" 4 

1.83970 x 10 -4 

1.46179 x 10~ 5 

1.336880384 x 10~ 4 

t = 

2.0 

7.1154 x 10" 4 

3.1006 x 10" 4 

1.35250 x 10 -4 

0.73805 x 10” 5 

9.818252567 x 10” 5 

t = 

3.0 

3.9213 x 10" 4 

1.7100 x 10 -4 

7.46014 x 10 -4 

0.40272 x 10' 5 

1.029870405 x 10~ 5 


Table 2 d: The rate of convergence for At = 0.01 and At = 0.0001 respectively 


At = 
N = 

0.01 

50 

Present (p = 1) 
3.9213 x 10" 4 

order 

At = 0.0001 
N = 50 

Present (p = 1) 
3.9213 x 10 -4 

order 

N = 

100 

9.9430 x 10" 5 

1.9773 

N = 100 

9.9430 x 10“ 5 

1.9909 

N = 

150 

4.4895 x 10" 5 

1.9649 

N = 150 

4.4895 x 10“ 5 

2.0032 

N = 

200 

2.5808 x 10" 5 

1.9192 

N = 200 

2.5808 x 10“ 5 

1.9931 

N = 

250 

1.6965 x 10" 5 

1.8734 

N = 250 

1.6965 x 10“ 5 

1.9929 


The corresponding graphical illustrations are presented in Figures 2 for k\ = —2, k 2 — 1, 
k 3 = 1, N = 400 and At = 0.001 at different t for best parameter p = 0.0002166. In 
Figure 3-4, computed solutions of v different time levels for k t ,k 2 fixed and ki,k 3 fixed 
respectively. 


Figure 2: Numerical Solutions at various 
t for N = 400, At = 0.001, p = 0.0002166 
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a: k\ = — 2 , k ‘2 = 1, k 3 — —8 b: k\ = — 2 , k 2 = 1, fc 3 = —4 c: = — 2 , /c 2 = 1, A ;,3 = —0 
Figure 3: Computed solutions of V Problem 1 for different time levels (ki, k 2 fixed) 


a: k 


1 = -2, k 2 = - 8 , k 3 = 1 b: ki = -2, k 2 = -4, k 3 = 1 c: £4 = -2, /c 2 = 0, k 3 
Figure 4: Computed solutions of C Problem 1 for different time levels (ki, k 3 fixed) 


1 


a: k\ = — 8 , k 2 — 1, k 3 = 1 b: = —4, k 2 — 1, = 1 c: hi =0, k 2 — 1, fc 3 = 1 

Figure 5: Computed solutions of V Problem 1 for different time levels (k 2 , k 3 fixed) 

Problem 2) Numerical solutions of considered coupled Burgers’ equations are obtained 
for k\ — 2 with different values of k 2 and k 3 at different time levels. In this situation the 
exact solution is 


U (x, t) — a 0 — 2A( 4 f 2 2 k3 ^ 1 ) tanh(^4(x - 2At)) 

V(x,t) = a 0 (H*Ei) - 2 ^( 4 fc 2 ^-i ) tanh ( A ( x ~ 2At )) 


10 



Thus, the initial and boundary conditions are taken from the exact solution is 


U(x, 0) = a 0 - 2A( 4 f 2 l^ 1 ) tanh(Ax) 


V ( x , 0) = a 0 (|g_i) _ tanh(^x) 


Thus, the initial and boundary conditions are exracted from the exact solution. Where 
ao = 0.05 and A = — (.The numerical solutions have been computed for the 


domain x G [0,1], At = 0.001 and number of partitions N = 10 and 100. The maximum 
absolute errors have been computed and compared in Tables 3 a-3 b for t — 1 with those 
available in the literature [25]. 


Table 3 a: Error norms for t — 1, At — 0.001, U(x,t), k\ = 2, k 2 = 1 and = 0.3 

Present (p = 1) [25], (A = 0) [25] (Various A) 

N = 10 3.7323 x 10~ 6 3.73505 x 10~ 5 0.00077 x 10“ 5 (A = 6 x 10" 5 ) 

N = 100 3.7350 x 10~ 6 3.73503 x 10" 5 0.00078 x 10" 5 (A = -4.087 x 10" 5 ) 


Table 3 b: Error norms for t — 1, At — 0.001, V(x,t), k\ = 2, k 2 = 1 and k% = 0.3 

Present (p = 1) ^5], (A = 0) [22] (Various A) 

N = 10 1.2569 x 10" 6 1.29030 x 10~ 5 0.00079 x 10" 5 (A = -6 x 10~ 5 ) 

N = 100 1.2871 x 10~ 6 1.29038 x 10~ 5 0.00079 x 10“ 5 (A = -4.087 x 10" 4 ) 


Table 3 c: 

Maximum error norms 

for U(x,t ) in Problem 2 (N = 21, A t = 0.01, k\ = 2) 

t 

k 2 

co 

Present (p — 1) 

pm 

ra 

0.5 

0.1 

0.3 

8.8160 x 10“ 6 

4.173 x 10~ 5 

4.189217417 x 10~ 5 


0.3 

0.03 

9.2556 x 10“ 6 

4.585 x 10~ 5 

4.584830094 x 10“ 5 

1.0 

0.1 

0.3 

8.8878 x 10" 6 

8.275 x 10" 5 

8.269641708 x 10“ 5 


0.3 

0.03 

9.3324 x 10" 6 

9.167 x 10" 5 

9.147335667 x 10” 5 

3.0 

0.1 

0.3 

8.9174 x 10" 6 

2.408 x 10” 4 

2.401202768 x 10" 4 


0.1 

0.03 

9.3691 x 10" 6 

2.747 x 10~ 4 

2.704203611 x 10" 4 


Table 3 d: 

Maximum error norms 

for V(x,t) in Problem 2 (N = 21, A t = 0.01, k\ = 2) 

t 

k 2 

h 

Present ip — 1) 

m 

ng 

0.5 

0.1 

0.3 

2.8380 x 10“ 6 

5.418 x 10~ 5 

9.094743099 x 10~ 6 


0.3 

0.03 

1.1179 x 10" 5 

2.826 x 10~ 5 

2.48218881 x 10" 5 

1.0 

0.1 

0.3 

2.8686 x 10“ 6 

1.074 x 10" 4 

1.696286567 x 10" 5 


0.3 

0.03 

1.1269 x 10" 5 

5.673 x 10" 5 

4.965329678 x 10" 5 

3.0 

0.1 

0.3 

2.9081 x 10" 6 

3.119 x 10" 4 

4.505480184 x 10“ 5 


0.1 

0.03 

1.1301 x 10" 5 

1.663 x 10" 4 

1.498311672 x 10“ 5 
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Figure 6: Numerical Solutions for U(x,t ) and V(x,t), N = 21, A t = 0.001, t — 1, 0 < x < 1, k 2 — 0.1 and 

Problem 3) Consider the Coupled Burger Equation system (|T|) with the following initial 
conditions 


f/(x,0) = | 

V{x,0) = l 


sin(27rx), 

0, 

0, 

— sin(27rx), 


x G [0,0.5] 
x G (0.5,1] 

x G [0,0.5] 
x G (0.5,1] 


and zero boundary conditions. In the Problem 3, the solutions have been carried out on 
x G [0,1] with At = 0.001 and number of partitions as 50. Maximum values of u and v 
at different time levels for k 2 = k 5 = 10 have been given in Table 4a and 4 b, while the 
Tables 4 c and 4 d represent the maximum values for k 2 = k% = 100. 


Table 4 a: Maximum values of U at different time levels for k 2 = k% = 10 


t 

Present (p — 1) 

m 

pa 

at point 

0.1 

0.144501 

0.14456 

0.144491495800 

0.58 

0.2 

0.052353 

0.05237 

0.052356151890 

0.54 

0.3 

0.019317 

0.01932 

0.019318838080 

0.52 

0.4 

0.007183 

0.00718 

0.007184856672 

0.50 


Table 4 b: Maximum values of V at different time levels for k 2 = k% = 10 


t 

Present (p — 1) 

m 

Pi 

at point 

0.1 

0.143155 

0.14306 

0.143141957500 

0.66 

0.2 

0.047003 

0.04697 

0.047006446750 

0.56 

0.3 

0.017258 

0.01725 

0.017260356430 

0.52 

0.4 

0.006415 

0.00641 

0.006416614856 

0.50 
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Table 4 c: Maximum values of U at different time levels for k 2 = k 3 = 100 


t 

Present (p = 1) 

m 

ra 

at point 

0.1 

0.04168 

0.04175 

0.041682987260 

0.46 

0.2 

0.01476 

0.01479 

0.014770415340 

0.58 

0.3 

0.00533 

0.00534 

0.005337325631 

0.54 

0.4 

0.00197 

0.00198 

0.001978065014 

0.52 


Table 4 d: Maximum values of V at different time levels for k 2 = k 3 = 100 


t 

Present (p — 1) 

m 

ra 

at point 

0.1 

0.05074 

0.05065 

0.050737669860 

0.76 

0.2 

0.01035 

0.01033 

0.010356602970 

0.64 

0.3 

0.00351 

0.00350 

0.003517189432 

0.56 

0.4 

0.00129 

0.00129 

0.001294450199 

0.52 


Figs. 7, 8 and 9 show the numerical results obtained for different time levels t G [0,1] 
at k 2 = k 3 = 10 for U and V with different values of k\. From the Figs. 7-9, it can be 
easily seen that the numerical solutions U n and V n decay to zero as t and k\ increased. 


Fig 7: Num. Sol. U(x,t) and V(x,t) of Problem 3 at different time levels for k 2 = k 3 = 10 while k\ — 1 


Fig 8: Num. Sol. U(x,t ) and V(x,t) of Problem 3 at different time levels for k 2 = k 3 = 10 while k\ = 10 
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Fig 9: Num. Sol. U(x,t) and V(x,t) of Problem 3 at different time levels for &2 = ks = 10 while k\ = 50 

5 Conclusion 

The collocation method together with the exponential B-spline as trial functions has pre¬ 
sented to get the numerical solutions of the coupled Burgers’ equation system. The free 
parameter in the exponential B-splines is searched experimentally to get the best numer¬ 
ical solution for the first problem. In the other problem, results are documented for p — 1 
as an example. The proposed method has produced less error than the methods listed in 
the tables for some text problems. The results are satisfactory and competent with some 
available solutions in the literature. Another advantages is that the method can be used 
without the complex calculations, to solve the system of differential equations reliably. 

And also the collocation method together with B-spline approximations represents an 
economical alternative since it only requires the evaluation of the unknown parameters at 
the grid points. 
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